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Abstrac t. A large-scale three-dime nsional model of Galactic extinction is presented based on the Galactic dust distribution 
model o f iDrimmel & Soergell JioOll) . The extinction Ay to any point within the Galactic disk can be quickly deduced using a 
set of three-dimensional cartesian grids. Extinctions from the model are compared to empirical extinction measures, including 
lines-of-sight in and near the Galactic plane using optical and NIR extinction measures; in particular we show how extinction 
can be derived from NIR color-magnitude diagrams in the Galactic plane to a distance of 8 kiloparsec. 
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1. Introduction 

In the past a fundamental obstacle to Galactic studies has 
been extinction due to interstellar dust, which has limited our 
view of the Galactic stellar distribution to the solar neigh- 
borhood. During the past decade NIR surveys, together with 
ever deeper optical surveys, have been piercing the inter- 
stellar haze, probing the stellar distribution in the Galactic 
plane and revealing its nonaxisymmetric structure, including 
the bar red nature of the Galactic bulge ( Blitz & Soerael 199l| 
fWeiland ltai]ll994l). the Gala ctic warp jPiorgovsk i & Sosin 
1989; Freude nreich et alJll994l) and the spiral arms ( Drimmel 
i:2000). Newly available high-resolution NIR surveys, such as 
2MASS and DENIS, promise to further reveal the nature of 
these structures. However, while the effects of interstellar dust 
are mitigated at NIR wavelengths, they are still important. 

The most recent map of Galactic extinction over the entire 
sky is that of Schlegel, Finkbeiner and Davis (1998, hereafter 
SFD), however it only maps the total Galactic extinction and is 
therefore most appropriate for extragalactic studies. To study 
stellar populations or objects found throughout the Galactic 
disk a volumetric description of Galactic extinction is required. 
Three-dimensional extinction models have been produced us- 
ing reddening data from stellar samples (Neckeletal. 1980; 
Arenou et al. 1992; Hakkila et al. 1997; Men de z & van Alteria 



1998h . but these give only a local description of extinction, 
being reliable to heliocentric distances of about 2 kiloparsecs 
due to limited sampling at optical wavelengths. However, in 
the NIR a well-defined stellar population can provide extinc- 
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tion measures along a line-of-sight to much greater distances; 
here we describe in detail one such method, deriving extinction 
measures from NIR color-magnitude diagrams (CMDs) using 
the known properties of the red-clump giants. This method is 
capable of rendering extinctions to distances as far as 8 kpc in 
the Galactic plane. 

In lieu of empirical extinction measures, a three- 
dimensional Galactic dust distribution model can be adopted 
from which extinction is derived. Previous studies at high and 
mid galactic latitudes have used "sla b" models, where a ver- 
tical dust density profile is adopted. IChen etal J (119991) have 
extended the usefulness and accuracy of such a model to rel- 
atively low latitudes by renormalizing their model to the SFD 
extinction map. Large-scale models that describe the dust dis- 
tribution over the entire Galactic disk have in the past been 
axisymmetric, coiTelating the radial variation of the dust dis- 
tribution with that of the gas (hydrogen) surface density, de- 
rived from HI, HII and CO observations together with as- 
sumed values an d gradients f or the gas:dus t (metallicity) an d 
CO:H2 ratios (eg. lWainscoatet al.. 1992; Ortiz & Lepinell993l) . 
Recently models have been constructed based upon FIR data, 
where Galactic emission is dominated by the thermally ra- 
diatin g dust. Initially these models wer e also axisymmetric 
( Sper gel et alJ Il996t iDavies et all Il997l) . but more recently 
have adopted nonaxisymmetric structures to account for their 
obser ved FIR emission, including the Galactic warp an d spiral 
arms iDrimmel & Spergei200lllLaunhardt et alJ2002l) . 

Here we present for general use a large-scale three- 
dimensional model of Galactic extinction based on the Galactic 
dust distribution model of .Drimmel & Spergel C200L hereafter 
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DSOl). The following two sections describe the construction 
and application of the extinction model. In section four we de- 
scribe how empirical extinction measures can be derived from 
NIR CMDs, and we compare the predictions of the model with 
these and other extinction measures. The summary section de- 
tails a number of caveats that a potential user of the model 
should be aware of. 

2. Galactic extinction model 

Given a three-dimensional model of the dust distribution, pd(x) 
the extinction in a given direction can be found by integrating 
along a line-of-sight to the point of interest, 

Av(x) = 1.086rv = 1.086/fv f Pd ds, (1) 

Jo 

while the total Galactic extinction is found by integrating to in- 
finity. The model of the dust distribution is detailed in DSOl. 
We mention here that this model is composed of three struc- 
tural components: a warped, but otherwise axisymmetric disk 
with a radial temperature gradient, spiral arms as mapped by 
known HII regions, and a local Orion-Cygnus arm segment. 
The structural parameters of the dust model are constrained by 
the FIR observations of the COBE/DIRBE instrument, while 
the determination of the parameter /cy, and the decomposition 
of the flux density into dust density and emissivity, is achieved 
by modeling the extinction in the COBE/DIRBE NIR observa- 
tions. However, the model as presented here has been amended 
in its description of the spiral arm geometry. 

The previous spiral arm geometry used in DSOl 
was based on the mapping of Galactic HII regions 
llGeorgehn & GeorgeUrJ[l97ft iTavlor & Cordeslll993h . Due to 
a lack of data, this map did not extend to the opposite side of the 
Galaxy. In order to provide a realistic model of the spiral arm 
contribu tion over the whole of the Galaxy, t he extended geom- 
etry of Blan d-Hawthorn & MalonevI (?2002') is adopted. Fig. ^ 
shows the dust surface density of the model, including the new 
spiral arm component. No other parameters of the dust model 
have been changed, and in fact the extension of the arms does 
not lead to any significant differences in the FIR emission pro- 
files as predicted by the model of DSO 1 , because the spiral arms 
on the far-side of the Galaxy are unresolved in the FIR DIRBE 
data. This model will need further refinement once FIR emis- 
sion toward the Galactic center is considered, as it predicts arm 
tangents within |/| < 20° which are not evident in the Galactic 
plane FIR emission profile. 

Another feature of the dust model of DSOl is the use of 
direction-dependent rescaling factors that are based on the FIR 
residuals between the DIRBE 240 yum data and the predicted 
emission of the parametric dust distribution model, and effec- 
tively "correct" the dust column density of the smooth model to 
account for small (angular) scale structure not described explic- 
itly by the model. In practice, for any given DIRBE pixel, one 
of the three structural components is rescaled to reproduce the 
FIR flux. The component chosen for rescaling is typically that 
which needs the least fractional change in its column density to 
account for the FIR residual, though the spiral arms are prefer- 
entially chosen near the Galactic plane (see DSOl for details). 




Fig. 1. Map of the dust surface density from the dust distribu- 
tion model with the extended geometry for the spiral arms. The 
Sun's location is indicated by the point at center- left. 

The effect of the rescaling factors is to add (detract) dust along 
the entire line-of-sight /or the rescaled component only, there- 
fore their application represents an approximate correction to 
the model in the sense that no bias with respect to distance is 
eff'ected, apart from that described by the rescaled component. 
It should also be noted that the rescaling procedure implicitly 
assumes that the deviation between the predicted and observed 
FIR emission is not due to variations in dust temperature. 

In DSOl the rescaling factors were used to refine the dust 
model when accounting for extinction in their modeling of the 
J and K band emission observed by COBE, hence only applied 
for galactic latitudes \b\ < 30°. As a practical matter, the rescal- 
ing factors for latitudes \b\ > 30° are based on the SFD Galactic 
extinction map, as the low signal-to-noise of the 240 /urn emis- 
sion at high galactic latitudes does not allow their reliable de- 
termination. Rescaling to the SFD extinction map was done as 
follows: The rescaled dust density can be expressed as 

Pd = ^ fiPi (2) 

where the sum is over the structural components of the dust 
density model (/ - disk, spiral arms, local Orion arm), and 
the rescaling factors /j are direction dependent, i.e. fi(l,b). 
Rescaling is applied to only one component, so that at most 
one of the factors /i 7^ 1 for a given line-of-sight. For any given 
direction the total extinction is then 

Av(cx,) = 2/;A,(tx)), (3) 
where 

Pi ds, (4) 
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the integral being taken along the line-of-sight, p, correspond- 
ing to the dust density associated with component / and /cy is 
the mean opacity. Rescaling to the SFD extinction map, Asfd. 
is effected by insisting that A(oo) = Asfd> which leads to 

^ ^SFD - 2/*disk^/(°°) 

/disk = -. — — ■ (5) 

Adisk(oo) 

Only the disk component need be rescaled, as at high Galactic 
latitudes it is only this component that contributes significantly 
to the dust column density. 

Fig. 121 shows the skymaps of total Galactic extinction, both 
with and without rescaling. Rescaling effectively adds angular 
detail that is not described by the parametric dust model. To 
give the reader an appreciation of the importance of the rescal- 
ing factors at low galactic latitudes, J band emission profiles 
are shown in Figs. |3]and|4]that both employ and neglect their 
use. 

3. Application 

An estimate of the extinction to any point within the Galaxy 
using the rescaling factors can be achieved by either integrat- 
ing the rescaled dust model, pd, along the line-of-sight to the 
point of interest (Eq.[2with pd),or by summing over the extinc- 
tion contributions of each component, determined separately 
via a Une-of-sight integration, and applying the rescaling fac- 
tor to the appropriate component. However, such integrations 
are time consuming, especially if many lines-of-sight need to 
be considered. To provide a more efficient means of finding the 
extinction within the Galactic disk we have constructed a set of 
three-dimensional cartesian grids of extinction in V, each cor- 
responding to a separate component of the dust model, as well 
as a table of the rescaling factors described above. Both large- 
scale and small scale grids are provided; the former cover the 
entire Galactic disk while the later describe the extinction about 
the Sun in greater detail. Using the three-dimensional extinc- 
tion grids a value for the extinction A, due to each of the com- 
ponents i of the dust distribution can be found for any point in 
the Galactic disk via interpolation. Together with the appropri- 
ate rescahng factor a final estimate of the extinction is arrived 
at: 

Ava/7,r) = 2 /;■(/, fe)A,(/,fe,r). (6) 

Trilinear interpolation from the grids described above is orders 
of magnitude faster than numerically integrating the dust model 
to individual sources, and the interpolation error is of the order 
of 0.2%. 

As an example. Fig. [S] shows a slice in the Galactic plane 
through the large-scale grid for the spiral arm component. Both 
this grid and that of the disk component have xyz dimensions 
of 30 X 30 X Ikpc, the z = plane corresponding to the Galactic 
plane. Thus they describe the extinction out to a galactocentric 
radius of R - 15kpc. The grid spacing parallel to the Galactic 
plane is 200 pc, while perpendicular to plane the resolution is 
20pc. Meanwhile, the large-scale grid for the Orion arm com- 
ponent has reduced dimensions and finer grid spacing parallel 
to the Galactic plane, so that it covers only part of the Galaxy, 
having dimensions of 3.75 x 7.5 x Ikpc. 




Fig. 5. Extinction in the Galactic plane due to the spiral compo- 
nent of the dust distribution out to R = 15kpc from the Galactic 
center. The Galactic center is at center, with the Sun located at 
center-left. The tangents to the spiral arms are clearly evident. 
The maximum extinction is 5.7 magnitudes in V. 



To provide greater detail near the Sun and to avoid inter- 
polation errors, two local grids are constructed, one for the 
Orion arm and another for the disk component; the spiral arm 
component does not make a significant contribution near the 
Sun. Both grids have their xy coordinates centered on the Sun, 
though with different grid spacing in the xy directions. (The 
grid spacing in the z direction for all grids is Az = 0.02 kpc.) 
The second of these local grids, that of the disk component, is 
necessary only to avoid interpolation errors which would other- 
wise be greater than 5% for heliocentric distances less than 500 
pc with the larger grid spacing of the large-scale grid. Ideally 
this region about the Sun should be described by a more de- 
tailed local model of the dust distribution. 

The rescaling factors over the entire sky are also sup- 
plied on the same resolution as the FIR COBE/DIRBE data: 
Because the rescaling factors are based on the residuals be- 
tween the DIRBE FIR observations and the predictions of 
the dust distribution model, the rescaling factors are rooted 
to the DIRBE data structure; for each DIRBE pixel there 
is an associated rescaling factor. This introduces some com- 
plication in the retrieval of the rescaling factors for any ar- 
bitrary direction as the DIRBE sky maps use a nonstan- 
dard projection and a binary pixel ordering scheme. (See 
http://space.gsfc.nasa.gov/astro/cobe/skymapinfo.html for a 
brief introduction and Appendix G of the COBE Guest 
Investigator Software (CGIS) Software User's Guide (version 
2.2), retrievable at ftp://rosette.gsfc.nasa.gov/pub/cobe-gi/doc/ 
for further details.) To apply the rescaling factors the user must 
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determine the DIRBE pixel that corresponds to the direction of 
interest. There are at least three options for this: 

1 . Write a routine that searches through the entire list of pixel 
coordinates for the nearest pixel. 

2. Work within the IDL environment developed by the 
COBE data reduction team which provides an ef- 
ficient means, via CGIS IDL code, to retrieve the 
pixel number nearest a given direction. This IDL pack- 
age with instructions for installation can be found at 
http://space.gsfc.nasa.gov/astro/cobe/cgis.html 

3. Use stand alone CGIS standard FORTRAN code to re- 
trieve the pixel number for a given direction. This 
code can be downloaded via anonymous ftp from 
'ftp://rosette.gsfc.nasa.gov/pub/cobe-gi/ 

Thus a file is provided for the user containing, for each DIRBE 
pixel, the DIRBE pixel number, its galactic coordinates, an in- 
dex of the component to be rescaled and its rescaling factor. 
This file, together with those containing the extinction grids 
and detailed instructions for their use, can be found at the 
anonymous ftp site ftp://ftp.to.astro.it/astrometria/extinction/ 

To summarize, an algorithm for finding the extinction Ay 
to a point in the Galaxy (/, b, r) is outlined below: 

- Find the COBE pixel coinciding to (I, b). 

- Recover the rescaling factor / and component / to be 
rescaled. (Set f - I for the other components.) 

- Determine the grid coordinates: 
{l,b,r)^{x,y,z)^{ij,k) 

- Interpolate from appropriate grids to find A,(/, b, r) for each 
component. 

- Sum over the components (Eq.|6} to arrive at Ay(l, b, r). 

An example of IDL and FORTRAN code that performs the above 
algorithm can be found on the anonymous ftp site mentioned 
above. 

It is important to note that if one decides not to use the 
rescaling factors it is sufficient to use only two grids, a large- 
scale and a local grid, which are each a sum of the grids men- 
tioned above, and perform only a single interpolation from the 
appropriate grid. Fig. |6l shows the resulting large scale grid, 
also provided to the user. One might wonder why the same is 
not done above, that is, construct a single three-dimensional 
grid of the final rescaled Av(/, j, k). The reason is that the use of 
the rescaUng factors would introduce discontinuities in such a 
grid, rendering interpolation unreliable; first interpolation must 
be done on each (smooth) component, then the direction de- 
pendent rescahng factor apphed. 

As a final detail we mention that the extinction from 
the model, given in the V band, can be transformed 
to other wavebands us ing the (A/Ay) ratios given by 
Rieke & Lebofskv ( 198^. Recently there has been some de- 
bate whetti^^ese ratios may be too high in the NIR ( Glass 
lT999l:lDrainel200,lh . however these ratios are preferred here 
in order to remain consistent with the NIR modeling of 
DSOl. It should also be noted that if using these ratios one 
assumes a specific reddening curve which in practic e is spa- 
tially variable blueward of the Johnson R band JMathid 
1X990; ■Fitzpatrick.1999.) . 
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Fig. 6. Extinction map in the Galactic plane resulting from all 
components of the dust model (disk, spiral arms and the local 
Orion-Cygnus arm). Orientation and scale is as in the previous 
diagram and the maximum extinction is 37 magnitudes in V. 

4. Comparing model with data 

The Galactic extinction map of SFD , Asfd, and of the extinc- 
tion model, Aoo, are both based in part on COBE FIR data, 
however the maps were derived using different approaches. It 
is therefore interesting to compare the two extinction maps 
before using the SFD map to rescale the high latitude data. 
The two maps of Galactic extinction were compared for each 
COBE pixel between galactic latitudes 5° < \b\ < 40°. The 
lo wer boundary i s the limit of reliability for Asfd as stated 
bv ISchlegel etalJ flQQS). while above \b\ = 40° the 240|/m 
data becomes too noisy to derive reliable rescaling factors. 
In addition data in the directions of the Orion nebula, Rho 
Ophiucus, Andromeda and the LMC and SMC were removed. 
The mean residual of the remaining data, (Asfd - Aoo), was 
found to be 0.096 magnitudes, while the mean relative differ- 
ence, ((Asfd ~ ^c»)/Asfd)i was equal to 0.024, showing that 
the two maps have nearly the same normalization. This con- 
cordance between the two extinction maps is the reason why 
the boundary for the different rescaling schemes at \b\ - 30° 
is hardly evident in Fig.|2] However, the mean absolute resid- 
ual, (|AsFD - AooD has a value of 0.15 magnitudes, or a mean 
relative absolute difference of 0.19, showing that there is sub- 
stantial scatter in the difference between the two maps. 

However, rather than provide a sky map of the to- 
tal Galactic extinction, the primary purpose of the extinc- 
tion model is to give an estimate of extinction to any 
point within the Galactic disk. To evaluate the model's 
performance in this regime we compare extinction esti- 
mates from the model with empirical extinction measures 
using NIR data from the second incremental data release 
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Fig. 7. NIR color-magnitude diagram for a typical low-latitude 
field (/ = 65°, b = 0°), taken from the 2MASS survey. Dashed 
curves delineate the selected region isolating the red-clump gi- 
ants. Filled circles show the maxima of the red-clump for indi- 
vidual magnitude bins. 
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Fig. 8. Gaussian fits (lines) to the red-clump counts (points) in 
three magnitude bins for the field / - 65° b - 0°. Solid line 
(filled circles) corresponds to the 12.7<mK<13 bin, the dashed 
line (pluses) to the 12.1<mK<12.4 bin, and the dot-dashed line 
(open circles) to the 11.5<mK<11.8 bin 



(2IDR) of the 2MASSi project ('Skrutskie et alJ ( ll997lEoO0l) . 
[http ://www.ipac.caltech.edu/2mass/releases/docs.htmll. 

In iLopez-Corredoira et al 1 (l2002i hereafter L02) a method 
was presented to obtain both the star density and interstellar ex- 
tinction along a line-of-sight by extracting the well-known red- 
clump population (spectral type K2III) from the infrared color- 
magnitude diagram (CMD). These stars constitute the majority 
of the disc giants (Cohen et al. 2000; Hammersley et al. 2000) 
and can be easily identified in the NIR CMDs. The method is 
extensively described in L02, so we give only a brief summary 
with some additional details here. 

(J - K, otk) CMDs are built for 1° x 1° fields by using avail- 
able 2MASS data. In these diagrams stars of the same spec- 
tral type (i.e. the same absolute magnitude) will be placed at 
different locations in the CMD; the effect of distance and ex- 
tinction cause the red-clump giants to form a broad diagonal 
branch running from top left to bottom right in the CMDs. In 
order to isolate the red-clump sources in the CMD we use the- 
oretical traces o f diff'erent spectral type s, based on the updated 
"SKY" model dWainscoat et alJll992l) . to define the possible 
color range of the K-giant branch in the CMDs, without any 
further implication in the method. (The "SKY" model uses a 
double exponential disk for the dust distribution, used here to 
roughly approximate interstellar extinction.) For each field ap- 
propriate traces are chosen to isolate the K-giants and avoid 
contamination by other stellar populations, especially dwarf 
stars and M-giants (see Fig.Q. 

Once the optimal traces have been selected, the giant stars 
are extracted from the CMD and binned in apparent K magni- 
tude, otk. For each magnitude bin, count histograms in color are 
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constructed. A Gaussian function was then fit to the histograms 
to determine the color of the peak counts at each magnitude: 

/(x;A,;U,cr) = Aexp[-(x-;u)V2o-2], (7) 

where x is the binned (J-K) color and /i is taken as the color of 
the peak. Fig.|8lshows the Gaussian fits to the color histograms 
obtained at three different magnitude bins for the field of Fig. 
The maxima are identified as corresponding to the mean ob- 
served color (J-K)„,^ of the K2III stars at the given apparent 
magnitude otk since they are by far the most prominent popu- 
lation in the selected region. (See also Fig. 2 in L02 for further 
details). 

The extinction Ak(hjk) can be determined by tracing how 
the peak (J-K) of the red-clump counts changes with otk and 
the intrinsic mean c olor (J-K)o of the red-clu mp. From the 
color excess and after lRieke & Lebofskvl ( ll985l) : 

(J-K)„,-(J-K)o 

Ak('«k) = ■ 

Once the extinction Ak is estimated at the apparent magnitude 
Otk, a mean distance to the stars can then be found given their 
mean absolute magnitude. For the K2III population the mean 
absolute magnitude and the intrinsic color were assumed to 
be Mk = -1.65 and (J-K)o=0.75, with a Gaussian disper- 
sion of 0.3 mag in absolute magnitude and 0.2 in color (see 
L02). Recent results in open clusters yield very similar values: 
Grochols ki & Sara jedini ( 2002. ) obtained Mk = -1.62 with a 
standard deviation of 0.21 and mean color (J-K)o of around 0.7 
for the old disc population, with small dispersion due to metal- 
licity or age gradients, while . Pie trzvnski et al. (2003) show that 
the K-band magnitude of the red-clump is a good distance indi- 
cator, independent of the metallicity and age, with some small 
dispersion in the color dependence. 
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Uncertainties of the extinction were estimated from the un- 
certainty of the mean color of the Gaussian component; 

(TA, =cr/(1.52VA^) (9) 

where cr is the standard deviation of the fitted Gaussian and 
A' — Act ^^ITjt is the estimated number of red-clump stars. This 
uncertainty gives an indication of the statistical error only, and 
does not include other possible systematic errors. Some uncer- 
tainties of the method, such as metallicity effects for the clump 
giants or differences between the real red-clump color distribu- 
tion and the assumed Gaussian function are fully described in 
L02, so they will not be repeated here. 

Critical aspects that are to be considered in the Gaussian fit 
are the possible contamination by dwarf stars and the effect of 
the completeness limit of the 2MASS survey. For fields such 
as those used here, dwarf contamination is negligible for mK ~ 
13 (see §3.3 in L02 for details) while the 2MASS survey is 
complete up to magnitudes of mK ~ 14. For this reason, the 
Gaussian fits have been obtained in each field only up to mK ~ 
13, well above the completeness limit of the survey and where 
the dwarf contamination is small. 

This method has been used for near-plane regions in both 
the outer and inner Galactic disc (L02; Lopez-Corredoira et al. 
2003, in preparation) with very satisfactory results. Some es- 
timations of interstellar extinction have also been obtained in 
this way for the star co unt predictions of the Besangon model 
llRobin&Crezel ll986') in the Galactic Plane and the results 
shown there are in good agreement between the observed ex- 
tinctio n in the K-band and the model predictions dPicaud et alJ 
The main strength of the method is that it is empirical 
and we can extract the red-clump giants from the CMDs as- 
suming only their mean absolute magnitude and intrinsic color 

Figs El and [TO] show the extinction along selected lines-of- 
sight in and near the Galactic plane as given by the extinction 
model and the NIR data. Also shown where available are ex- 
tinction measures from individual OB stars from the Neckel 
et al. compilation (.Neckel et al. 1980. 1995). These extinctions 
are converted to the K ban d using (AkMv) = 0.1 12 as given by 
iRieke & Lebofskv ^('l985'). When present these give a measure 
of the extinction out to a couple of kiloparsecs from the Sun, 
whereas the NIR red-clump data give a measure of extinction 
as far as 8 kpc. In general the NIR extinction measures are more 
reliable as each measure (distance and extinction) is based on 
many stars, as described above. In addition the distances of the 
OB stars are spectro-photometric, and thus suffer from large 
uncertainties due to misclassification and cosmic variation. 

In Fig.|9llines-of-sight with and without OB stars are pre- 
sented; in those fields containing no OB stars (left column) 
there is good agreement between the model and the NIR ex- 
tinction measures, while in those fields containing many OB 
stars the model is consistently giving higher extinction than in- 
dicated by the NIR data. This can partly be understood by the 
effect of the OB stars on their environment, as they will heat the 
ambient dust leading to higher dust emission than is typically 
seen in the arms. This contributes to the erroneous rescaling 
seen in the field at Z = 173°. Also worth noting for this field, 
as well as that at / = 187°, is that the two data sets give incon- 
sistent measures of the extinction, which must monotonically 



increase with distance. This inconsistency might in part be ac- 
counted for if the extinction curve for these lines-of-sight were 
nonstandard, i.e. R\ - 5 rather than the canonical value of 3. 1 . 
In this case the extinctions of the OB stars would systemat- 
ically decrease by approximately 3/5ths while their distances 
would increase, giving much better agreement with the NIR 
data which is insensitive to variations in the extinction curve. 
In any case, the model falls between the two empirical mea- 
sures for these two fields. 

For the inner disk (Fig.llO> the deviations between the data 
and the model are not as easily understood. Typically these de- 
viations are of the order of 20%, though seem to be higher in 
directions corresponding to tangents of spiral arms (I - 50° and 
30°). Also shown for some of these lines-of-sight is the extinc- 
tion curve for a model in which the disk was rescaled rather 
than the nominally selected spiral arms. The performance of 
this alternative rescaling is in general not significantly worse 
nor better, though it seems to be better along spiral arm tangents 
and worse in other directions. The exception is in the direction 
I - 50°, which is much better. However, this apparent improve- 
ment may be misleading, as the failure of the model may lie in 
its assumed geometry of the Sag-Car spiral arm, which places 
its tangent inside of I - 50°. The direction / = 75° is influ- 
enced by the local Orion-Cygnus arm, but shows good agree- 
ment with the NIR data out to a distance of 4 kpc, at which 
point the NIR data shows a large step in extinction. (The non- 
scaled extinctions suffer artifacts from interpolation errors as 
the single grid containing the dust model does not have suffi- 
ciently small gridding to resolve this feature.) Again, for this 
field and that at I - 60°, the OB extinctions are higher than 
those given by the NIR data. 

5. Summary 

We have here presented a three-dimensional Galactic extinction 
model, based on the dust distribution model of DSOl, which 
may be used for a variety of studies involving distant targets 
located in the Galactic disk or beyond. See Belokurov & Evan^ 
(I2OO3I) as an example. 

Any potential user of the extinction model should be aware 
of a number of caveats: 

1 . Absorptions to points within a few hundred parsecs of the 
Sun will not be reliable as the spatial resolution of the dust 
model does not allow a detailed description of dust in the 
vicinity of the Sun. A more detailed local map of extinction 
is needed for nearby sources. 

2. The dust model, even with rescaling, has finite angular res- 
olution. The size of the COBE pixels are approximately 
0.35° X 0.35°, with some variation over the sky. However, 
any field of this size (or even smaller) will contain unre- 
solved structures in the dust distribution. Thus the run of 
extinction as given by the model for a given line-of-sight 
should be considered as a mean extinction. 

3. The dust model does not include features associated with 
the Galactic bar nor the nuclear disk, though a central 
"hole" in the dust distribution is included. However, the 
parameters describing this hole are ad hoc as the FIR data 
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within 20° of the Galactic center was not modeled by DSO 1 . 
In addition, the adopted spiral arm geometry extend the 
arms too far into the center of the Galaxy, leading to arm 
tangents that are not evident in FIR emission in the Galactic 
plane. The extinctions within 0.35Rq of the Galactic center 
may suffer significant systematic errors. 

4. The density of the spiral arms are not well constrained by 
the COBE NIR data because the arm tangents are too dis- 
tant to cause strong extinction features in the diffuse NIR 
radiation; lines-of-sight corresponding to arm tangents may 
have large systematic errors. 

5. In directions where extragalactic sources contribute signif- 
icantly to the FIR flux the derived rescaling factors will be 
erroneous. This includes the Andromeda Galaxy, M33, the 
LMC and SMC. 

6. The rescaling factors in directions which contain anoma- 
lous dust temperatures will be in error. These include the 
nearby star forming regions of the Orion and rho Ophiucus 
complexes. 

For the purpose of validating the extinction model we 
have also presented a method by which the extinction can be 
empirically estimated using NIR color-magnitude diagrams. 
Comparisons made for lines-of-sight in the outer Galaxy show 
good agreement between the model and the NIR data, show- 
ing differences less than 0.05 magnitude, though directions to- 
ward known OB associations have a high relative difference, 
in part due to anamolous dust temperatures. Toward the in- 
ner Galaxy there is less concordance, especially in directions 
corresponding to spiral arm tangents. Most of our compar- 
isons with the NIR data are in the Galactic plane {b = 0). 
Interestingly the effect of the rescaling factors in the Galactic 
plane is not as clearly advantageous as it is at low Galactic 
latitudes {\b\ » 2 - 3°). Again, anomolous dust temperatures 
associated with star forming regions may be Umiting the use- 
fuUness of the rescaling factors. 

When reliable extinction data is available it is always 
preferable to rely on these rather than a model. However, even 
when such data is available a model may be useful to extrapo- 
late beyond the limits of the data, as it may be able to provide 
extinctions at larger (or smaller) distances or at higher resolu- 
tion than the data. The user is also free to compute new rescal- 
ing factors, or to rescale the entire model, when it is in clear 
conflict with data for a particular line-of-sight. 

Some of the caveats listed above suggest possible improve- 
ments that could be made in the Galactic dust distribution 
model, particularly toward the Galactic center. Work is cur- 
rently underway toward the development of a dust distribution 
model that also models the FIR emission within 20° of the 
Galactic center, not originally considered in the modeling of 
DSOl . With the recent availabiUty of high quality NIR data over 
the whole sky, together with methods for estimating the extinc- 
tion along lines-of-sight close to the Galactic plane, it will soon 
be possible to improve the extinction model by providing fur- 
ther observational constraints. Such added constraints will be 
particularly helpful in deriving a more reliable estimate of the 
contribution of the spiral arms. 
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Fig. 2. Sky maps of the predicted total Galactic extinction in a MoUweide projection, without (top) and with (bottom) rescaling. 
The direction to the Galactic center is at the center of the plot, positive galactic longitudes to the left. The quantity plotted in grey 
scale is - log(Av). 
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Fig. 3. J band emission profiles at low galactic latitudes, as observed by the COBE satelUte (crosses) and predicted by the dust 
model using rescaling factors (from DSOl). 
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. 4. As in previous figure, but using the dust model without rescahng factors. 
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Fig. 9. K band extinctions along lines-of-sight in the outer disk of the Galaxy. Solid and dashed lines show the extinction 
according to the model presented here, with and without rescaling factors respectively. Diamonds show extinctions based on NIR 
CMDs of 1° X 1° fields using 2MASS data. Asterisks show extinction measures from individual OB stars (Ne ckel et aL.199^ . 
the large asterisks corresponding to stars within 0.5° of the directions indicated, small asterisks stars within 1°. 
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Fig. 10. K band extinctions along lines-of-sight in the inner disk of the Galaxy. Lines and symbols as in previous figure, with 
the addition of a dotted line showing the extinction of the model in the case where the disk component of the model is rescaled 
rather than the spiral arm component. 



